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A large area (128m 2 ) Muon Tracking Detector (MTD), located within the KASCADE experiment, has been 
built with the aim to identify muons (S„ >0.8GeV) and their directions in extensive air showers by track 
measurements under more than 18 r.l. shielding. The orientation of the muon track with respect to the 
shower axis is expressed in terms of the radial- and tangential angles. By means of triangulation the muon 
production height is determined. By means of Hp, a transition from light to heavy cosmic ray primary 
particle with increasing shower energy E from 1-10 PeV is observed. Muon pseudorapidity distributions for 
the first interactions above 15km are studied and compared to Monte Carlo simulations. 



1. Introduction 

Muons have never been used up to now to recon- 
struct the hadron longitudinal development of EAS 
with sufficient accuracy, due to the difficulty of build- 
ing large area ground-based muon telescopes Q]. 
Muons are produced mainly by the decay of charged 
pions and kaons in a wide energy range. They are not 
always produced directly on the shower axis. Multi- 
ple Coulomb scattering in the atmosphere and in the 
detector shielding may change the muon direction. It 
is evident that the reconstruction of the longitudinal 
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development of the muon component by means of tri- 
angulation (2j E] provides a powerful tool for primary 
mass measurement [3. , giving an information similar 
to that obtained with the Fly's Eye experiment, but in 
the energy range not accessible by the detection of flu- 
orescence light. Muon tracking allows also the study 
of hadron interactions by means of the muon pseudo- 
rapidity [5]- Already in the past, analytical tools have 
been developed which describe the transformation be- 
tween shower observables recorded on the ground and 
observables which represent directly the longitudinal 
shower development [7] . Fig. 1 in ref [5] shows the ex- 
perimental environment. Measured core position dis- 
tributions for showers inside KASCADE range from 
40 - 120 m and inside Grande from 250 - 360 m.The 
shower core position ranges cover full trigger efficiency 
as confirmed by investigations of muon lateral density 
distributions [B] . With CR studies very high energies 
are accessible in the 'knee' energy region 10 15 - 10 16 5 
eV, which correspond to CM energies in the nucleon- 
nucleon system from 1.4-8 TeV currently covered by 
the Tevatron and the LHC. Provided, that it is feasible 
to focus on the first encounters with the atmospheric 
nuclei, muon multiplicity studies may provide insight 
into the high energy interaction at around 8 TeV. In 
this energy range saturation physics is expected to en- 
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Figure 1: p angle distributions (normalized to integral 
yield equal to one) for different muon size bins and 
different Ig(N^) /lg(N e ) ratio above 0.83 ('heavy') and 
below 0.83 ('light'), for 80 - 120 m core distance range 
and shower direction = 0° — 18°. Lines connect points. 



ter (rise of < px >) which is transformed to geometric 
scaling [2] by the color-glass-condensate theory. 



2. Muon Production Height 

The angular correlation of the muon tracks with 
respect to the shower axis is expressed by the p and 
the r angles pQ . The p angle contains some scattering 
which is represented by the r angle value exhibiting a 
a T ~ 0.2°. 

Fig. 1 shows p angle distributions for specific 
muon number lg{N^) bins corresponding to differ- 
ent shower energy bins @]. The p angle distribu- 
tions are plotted for 'light' and 'heavy' primary CR 
mass enriched showers, employing the lg(N lx )/lg(N e ) 
ratio (corrected for attenuation) [5] [3] to be larger 
('heavy') or smaller ('light') than 0.83. The distribu- 
tions show a dependence on the primary mass range, 
however, masked by the energy dependent penetra- 
tion. 



Based on p and t angles and the distance of the 
muon hit to the shower core i? M , the muon production 
height along the shower axis is calculated: 

hft = i? Al /tan(p- | r |) (1) 

The MTD-KASCADE system with its dense ar- 
ray grid and 80 - 120 m core distance range for the 
muon track allows to study the muon momenta in 
the 100 — 200 GeV range, hu will be considered 
for p > t which we can extend up to 20 km for a 
shower core muon hit distance window 80 - 120 m 
(note 100 m/20000 m = 0.28°) employing, accord- 
ing to simulations, ~ 200 GeV muons from above 
15 km. Fig. 2 shows muon production height dis- 
tributions for different muon size bins and different 
lg(N fJi )/lg(N e ) ratio above 0.83 ('heavy') and below 
0.83 ('light'). CORSIKA [12 simulations based on 
QGSjetII+FLUKA2002.4 (slope -2.7 and -3.1 be- 
low and above the knee, respectively) for Hydro- 
gen and Iron are shown in the Fig. 2 as well. In 
the low range the low energy interaction model 
(FLUKA2002.4) seems capable to describe the dis- 
tribution. The experimental distributions are getting 
more narrow with increasing energy but differently for 
'light' and 'heavy' CR primaries. The Snn numbers 
quote the CM energies assuming A=l CR primary. 
In the course of the analysis in |4] , /i M is transformed 
to [g/cm 2 ] and after subtracting from each 
an 'energy' dependent penetration depth the remain- 
ing depth exhibits [4] the mass A sensitivity. 
Mass resolution is limited and < >^< Hf[ > 
- 30g/cm 2 lnA [4]. 

Fig. 2 shows strong reduction of muons from 
above 15 km for 'light' primary CR particles with re- 
spect to Monte-Carlo especially for the highest en- 
ergy interval, which corresponds to a CM energy of 
~ 8 TeV. Standard Monte-Carlo (see above) predicts 
for a nucleon-nucleon collision charged particle mul- 
tiplicities < N ch > ~ 40 (2TeV) - 70 (8TeV). As a 
consequence, a fraction out of ~ 70 x 200GeV parent 
pions are missing. Very recent results from LHC at 
7 TeV report also [TU] smaller < N c h > at midrapi- 
tity compared to QGSjetll. It would be worthwhile to 
study to which extent the simulated distributions in 
Fig. 2 above ~ 15 km scale with the inelastic cross 
sections of H and Fe primary masses. At 10 16 eV on 
Nitrogen af nel amounts to ~ 400 mb and crf^ei ^ s about 
5 times larger. The ratio of charged particle multiplic- 
ities < N[ h e > / < N^ h > amounts to about 12. Both 
factors lead to the difference between H and Fe seen 
in muons in Fig. 2 at high altitude. Therefore, we as- 
sume that about one tenth out of 70 pions lead to an 
invisible energy fraction corresponding to 7 x 200GeV 
(~ ITeV) of the incident energy of - 10 4 5 TeV. For 
lower production height a regular shower development 
is taking over and described by the low energy inter- 
action model FLUKA2002.4. The slower development 
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Figure 2: Muon production height distributions 
(normalized to integral yield equal to one) for different 
muon size bins and different lg(N^)/lg{N e ) ratio above 
0.83 ('heavy') and below 0.83 ('light') and 80 - f 20 m 
core distance range and for shower direction 
6 = 0° - 18°. Also CORS1KA simulations based 
on QGSjetll + FLUKA2002.4 are given for Hydrogen 
(H) and Iron (Fe) CR primaries. Lines connect points. 



of showers seen in the 'light' data compared to simu- 
lated H, points in a classical picture to more ir decays 
in the simulations or more /i decays in the data, how- 
ever, can not be reached by decreasing Ui ne i in MC in 
the range as discussed in [TT] . 

The deviation from standard high energy Monte- 
Carlo is strongest for the 'light' CR primaries for 
which the quoted VSnn is around 8 TeV. The ob- 
servation that 'heavy' primaries show less deviation 
from the Fe prediction, points to a possible threshold 
effect for the 'light' CR primaries. 



3. Muon Pseudorapidity 

To investigate the deviation from the QGSjetll sim- 
ulations in Fig. 2 above 15 km especially for the high- 
est energy bin, pseudorapidity spectra for muon pro- 
duction height larger than 15 km have been cal- 
culated based on the p and r angles as introduced in 
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Figure 3: Pseudorapidity spectra (normalized to integral 
yield equal to one) for muon production height larger 
than 15 km and 80 - 120 m core distance range and for 
shower direction 6 = 0° - 18°. CORSIKA simulations 
based on QGSjetll + FLUKA2002.4 are given for 
Hydrogen (H) and Iron (Fe) CR primaries. Lines connect 
points. 

0- 

r, ~ ln(2p\\/pr) - -ln(Vp 2 +r 2 /2) (2) 

The pseudorapidity distributions are taken in the 
same size bins. The little variation of the n distri- 
butions is to some extent due to our analysis window 
(15 - 20 km). With our analysis window we filter out 
pions around their first interactions. Only very few 
'prompt' muons will be present from the very first in- 
teraction. We have to keep in mind that after a defi- 
nite number of interactions pion's probability of decay 
exceeds that of arriving at the next interaction level. 
One should note, that there is only little difference 
between 7/(pion) and ry(muon) [5 J. From independent 
studies with accelerators the effective beam energy for 
charged particle production (inelasticity K ~ 0.3) is 
quoted to be about 0.3 x Eq [13]. Effective beam ra- 
pidities %eam, are quoted in the Fig. 3 assuming A=l 
CR primary mass interacting in the atmosphere. 
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Good agreement for Hydrogen below the 'knee' and 
for Iron above the 'knee' is observed. Comparison of 
7] distribution for 'light' CR primaries and simulated 
H in the highest energy bin suffers from limited statis- 
tics. 

For the understanding of the r\ vea k positions far 
away from % eam the following considerations may ap- 
ply. When taking into account that the CR mass in- 
creases up to mass A ~ 56, beam rapidities are modi- 
fied by —0.5lnA, leading to very similar f)i> eam values. 
While the total energy increases from 10 15 — 10 165 
eV we know [3] [TS] (in the frame of QGSjet) that 
the mean mass of the CR primaries increases from 
< A >^ A to < A >^ 56. Therefore, correcting the 
quoted n beam gives the r) beam values. 

However these fjb ea m values exhibit a further shift 
to the r] pea k values possibly due to geometric scaling 

EH- 

Vpeak = (V(l + A)) x {f, beam - InA 1 ^), A ~ 0.2 (3) 

The n distributions are described by Monte-Carlo. 
Therefore, Monte-Carlo are able to cover the p\\ and 
Pt for the muons which stem mostly from charged pi- 
ons. Pions are considered to dominate at midrapidity, 
but here we deal with the pions in the fragmentation 
region which deliver muons conserving the rapidity 
of the parent mesons. With respect to the effect of 
missing muons at the highest energy, representing a 
separate feature, the comparison with other type of 
high energy interaction model would be of interest. If 
this effect is real, it should occur again at ~ 10 18 eV 
for 'heavy' CR particles. 

4. Conclusions 

Muon tracking allows to investigate ft,„ and r\ . Fu- 
ture analysis of other shower angle bins and a larger 
and improved quality data sample will provide a more 
detailed information on the nature of high energy 
shower muons. Also muon multiplicities provide valu- 
able parameters to support the study of the relative 
contributions of different primary cosmic ray parti- 
cles. A natural extension towards even larger shower 
energies is provided by KASCADE-Grande [15] . 
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